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Some of the new experimental results and theoretical developments presented at 
the Workshop on Deep Inelastic Scattering and Related Phenomena are reviewed. 



1 Introduction 

The present series of Workshops on deep inelastic scattering and related high- 
energy processes began in Durham in 1993. At that time results from the 
HERA ep collider were beginning to appear, and a forum where experimental- 
ists and theorists could get together to discuss these and other related measure- 
ments seemed appropriate. The meeting was a resounding success. The quality 
and quantity of the physics, together with the enthusiasm of the participants, 
all pointed towards the establishment of a 'deep inelastic scattering' workshop 
as an annual event. The Eilat (1994) and Paris (1995) meetings confirmed the 
Workshop as a truly international meeting, and one of the most important in 
the high energy physics calendar. This tradition has been continued in Rome, 
with a record number of participants and a wealth of interesting physics. 

Although HERA provided the original motivation, one of the keys to the 
success of the DIS Workshops is the way they draw together the whole deep in- 
clastic scattering community, with fixed-target experiments playing an equally 
important role. The hadron collider community has also been well represented, 
illustrating the complementarity of lepton-hadron and hadron- hadron collisions 
in providing information on hadron structure. 

There have been four main strands to the physics discussed at this Work- 
shop: (i) investigating the parton structure of the proton and photon as re- 
vealed in high-energy lepton-hadron and photon-hadron collisions respectively; 
(ii) understanding the origin of those events which are both deeply inelastic 
and diffractive; (iii) studying detailed QCD dynamics by means of particular 
hadronic final states (jets, heavy flavours, . . . ); and (iv) unravelling the spin 
structure of the nucleon by means of polarized deep inelastic scattering exper- 
iments. In this brief review I will attempt to highlight some of the new results 
in these different areas, together with their theoretical implications. The choice 
is necessarily restricted (lack of space and personal expertise being the main 
constraints) and does not come close to doing justice to all the interesting 
physics which has been presented and discussed. Nevertheless, I hope it will 
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give a flavour of what is without question one of the most important frontiers 
of particle physics today. 

2 Proton structure 

The traditional method of obtaining information on the parton structure of the 
nucleon is through measurements of deep inelastic structure functions. The key 
variables (for example, for the process e~ (fe) + p(p) — > e~ (k' — k — q) + X) are 
x and Q 2 , where 

Q 2 = -q 2 = -(k-k') 2 , x = Q 2 /2p-q. (1) 

Structure functions Ti(x, Q 2 ) are then obtained from the differential scattering 
cross section, d 2 a/dxdQ 2 . The main advance in recent years has been the 
dramatic increase in the range of x and Q 2 covered by experiment. With 
improvements in luminosity and detectors, HERA has been able to measure 
down to x ~ 10 -5 . At the same time, the Q 2 range has been extended 
both at the upper and lower ends. For t he former, this leads to an increase 
in quark substructure limits (see Section [^1]) and tests of perturbative QCD 
(pQCD) evolution. Data at low Q 2 are important for providing a bridge to the 
fixed-target data, and also for understanding the perturbative/nonperturbative 
transition. The region covered by the most recent HERA and fixed-target data 
is illustrated in Fig. 

2. 1 New structure function measurements 

Two significant improvements have been reported at this Workshop 0@. The 
1994 HERA F2 data now overlap with the fixed-target data, in the region 
Q 2 ~ 10 — 100 GeV 2 , x ~ 0.01 — 0.1, and the total ranges now spanned are 

1.5 < Q 2 < 5000 GeV 2 , 3 x 10" 5 < x < 0.3 . (2) 

As the kinematic region of the HERA F2 data continues to grow, two notable 
features persist: 

• F2 rises at small x for all Q 2 ; 

• NLO DGLAP evolution provides an excellent description of the Q 2 de- 
pendence (see next Section). 

The kinematic range of the ZEUS data has been further extended by the 
installation in 1995 of a beam pipe calorimeter. This allows electron and 
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Figure 1: x, Q phase space covered by various recent DIS experiments. 



positron scattering at smaller angles to be measured, which in turn leads, via 
Eq. (Q), to^smaller Q 2 values. Preliminary F-i measurements, reported at the 



Workshopu, are shown in Fig. Q The corresponding . 

ZEUS BPC: 0.16 < Q 2 < 0.57 GeV 2 , 



range is 
10~ 5 . 



(3) 



Fig. U also illustrates the complementarity of the HERA and fixed-target struc- 
ture- finction measurements. At the same low Q 2 values the E665 collabora- 
tions have measured for x ~ 10~ 3 — 10 -2 , see Fig. ||. All these data are 
important for constraining models of structure functions at low Q 2 . For ex- 
ample, the Donnachie-Landshoff Regge-based model B appears to give a good 
description of the Q 2 < 1 GeV 2 data, interpolating between the ZEUS and 
E665 data with a slowly rising ~ x~ 08 form. However this model does 
not include pQCD Q 2 evolution, and therefore fails to describe the steepening 
of F2 with increasing Q 2 which is apparent in the SVX/ISR data in Fig. |[ A 
comprehensive and critical review of models for F2 at low Q 2 can be found in 
Ref. □ Other new structure function data reported at the Workshop include up- 
dated measurements of F^ and Fg j F^ from the NMC collaborationEHl. These 
are particularly useful for constraining the medium- 2; quark distributions (in 
particular the u/d ratio) and will be incorporated in forthcoming 'global fit' 
analyses. 
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Figure 2: Preliminary ZEUS 1995 BPC F 2 measurements (top), together with ZEUS 1994 
SVX/ISR data at higher Q 2 (bottom), from-Ref. 0. The predictions of the theoretical models 

from Rcfs.rjQ are also shown. 



2.2 DGLAP evolution 



One of the central tenets of perturbative QCD is that the Q 2 evolution of 
structure functions J-i(x, Q 2 ) is determined by the DGLAP equations, provided 
that Q 2 is sufficiently large such that higher-twist (ex 1/Q 2 ) contributions 
can be neglected. More precisely, the theory predicts the factorization scale 
dependence of the quark and gluon distributions, q(x,/j, 2 ) and g(x,/j, 2 ), 

J_ / q \ oMj 2 } ( Pqq Pqg \ ( q \ 
" dn 2 [ 9 J 2n { P aq P a J [g ) ' W 
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Figure 3: F% measurements by the E665 collalmmjion @, with various theoretical model 

predictions lIlIu . 



where <g> denotes a convolution integral and the splitting functions have a 
perturbative expansion 

P ab (x, a.) = a s P^(x) + a 2 s P%>(x) + ... . (5) 

The structure functions are obtained as linear combinations of the parton dis- 
tributions: 

^(x,Q 2 ) = [ dy/y Y / q(y,Q 2 )®C q , i (x/y,a s (Q 2 )) 

+g(y,Q 2 )(g>C g 4x/ yi a s (Q 2 ))] . (6) 

For example, for Ti = x^F^ the coefficient functions are C q .i(z, a s ) — e 2 5(l — 
z) + 0(a s ), C g s(z,a s ) = 0(a s ). Once the distributions fi{x,Q 2 ) (i — q,g) 
and a s are specified |] at some 'starting' scale Qq, the theory predicts the 

a In practice the value of a s (M^) rather than a a (Q^) is used to quantify the strong 
coupling. 
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Figure 4: The Q 2 dependence of i<2 at small and medium i*', together with a NLO QCD fit, 

from the HI collaboration El . 



distributions at any Q 2 for which perturbation theory is valid. In practice the 
starting distributions are determined from the data by means of a global fit, see 
for example Refs. It is remarkable that only relatively simple functional 

forms are required, for example 



xfi{x, Ql) = A lX ~ x *(l + EiVS + - xY* 



(7) 



The splitting and coefficient functionsLare known exactly at leading and next- 
to-leading order (see for example Ref.Ellfor a compilation and list of references). 
Truncating at this order defines the NLO DGLAP system of equations. 

Several examples of NLO DGLAP fits to structure function data have 
been presented at the Workshop. Considering the wide range of processes and 
energies, the quality of the description is excellent. An example from the HI 
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collaboration L3 is shown in Fig. [|. What information can be extracted from 
such fits? At large x the DGLAP equation for F% reduces to dF^jd In Q 2 ~ 
ot s F2 ® Pqq, and a precision measurement of a s can be made, see Section || 
below. At small x the gluon is the dominant parton and so di^/^lnQ 2 ~ 
a s g £g) P qg - Thus while F2 directly measures the quark distribution at small x, 
its Q 2 variation measures the gluon. 

The shape of the starting distributions (Eq. (Q)) at small x is determined 
by the parameters Aj. The dominant partons here are the gluons and sea 
quarks, for which we may write 

xg{x, Ql) ~ A gX - x « , xq s {x,Ql) ~ A s x- Xs . (8) 

It is interesting to see how A g and As, obtained from fits to data, have evolved 
with time. Table |l| shows these parameters for various recent MRS anal- 
yses. |] Before HERA data became available, it was traditional to assume 
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Table 1: Evolution in time of the parameters determining the small-rr behaviour of the MRS 

quark and gluon distributions. 

Regge-motivated flat starting distributions (A « 0). Subsequently, theoret- 
ical studies of the BFKL equation for resumming large logarithms of 1/x 
(see Section |2.4| ) suggested that the behaviour could in fact be much steeper 
(A ~ 12a s \u2/tt ps 0.5). The early HERA data showed a somewhat less steep 
rise at small x, with fits giving A w 0.3. More recent HERA data show a pref- 
erence for different A5 and X g values. In fact in the most recent MRS fitsE3, 
where the minimum Q 2 of the fitted data is extended down to 1.5 GeV 2 , the 
gluon is 'valence-like' at the starting scale Q\ — 1 GeV 2 . Even more significant 
is the fact that the slope of the sea quark distribution is close to the Regge 

b Since all the fits listed have been performed at NLO in the MS scheme, the parameters 
corresponding to the same Qq are directly comparable. 
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prediction obtained from fits to the energy dependence of total hadronic cross 
sections B: 

xq S ~ a;l-a7»(0)«-0.08 at Q2 = 0(1 Gey2) _ (g) 

However the agreement between the fitted value of As and a-p(0) — 1 should 
not be taken too literally. The former is an unphysical parameter, depending 
to some extent on the factorization scheme, the value chosen for Qq, and the 
parametric form chosen for the starting distributions. Another problem with 
the Regge interpretation is that one would naively expect the gluon distribution 
to exhibit the same asymptotic behaviour for x — > 0. One can impose by hand 
the constraint Xs — A g , in which case the common value 0.04 (at Q% = 1 GeV 2 ) 
is obtainedE 2 ], but the quality of the fit is worse than when the parameters are 
allowed to be different. 

As Table [l] implies, the continual improvement in the structure function 
data requires a concomitant updating of the parton distributions derived from 
global fits. Thus 'new' data is able to discriminate between 'old' sets of partons. 
This is illustrated in Fig. |j| which shows recent HERA F2 data compared with 
the 1995 MRS £3 and 1QQ.4 GRVB 'dynamical parton' predictions. Also shown 
is the new MRS(R 2 ) fitl!3 which includes these data. Notice how the MRS(G) 
and GRV(94) predictions are now disfavoured - they predict too strong a Q 2 
dependence at small x. It will be interesting to see if the agreement between 
the dynamical parton model and the data can be restored by adjusting the 
parameters of the model, for example the value of the starting scale /z 2 , at 
which the distributions assume a valence-like form. 

We can summarize the results of this section by the statement that the 
small- a; structure function data are consistent with DGLAP evolution starting 
from a soft input for both quarks and gluons at a scale of order 1 GeV 2 . The 
fact that the input is soft means that we can obtain a simple analytic approx- 
imation to the solutions of the evolution equation in the (double) asymptotic 
limit Q 2 , 1/ x — > ooEj. In its simplest form, this result follows from the leading 
x — > behaviour of the lowest-order splitting function matrix: 

2 d ( q } a s (fi 2 ) ( \ ^ f q 



^{ 9 )-^r{2C F /x 2C A ,*)*[l)> ^ 
which leads, for sufficiently soft starting distributions, to the prediction0 

(11) 



exp 



C A a s , Q 2 ln Xq 



Ball and Forte have extended this 'double asymntntip scaling' result to include 
subleading corrections and NLO contributions O'ELS. An example of a com- 
parison between theory and experiment is shown in Fig. |(| Here F2 has been 
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Figure 5: The description of the rm*sti recent HERA data by &m 'old' parton distributions 
of Refs.HO, and the new MRS(R 2 ) fitlB. 



rescaled by a factor which is essentially the right-hand side of (|il|), computed 
at NLO, in order to remove the leading asymptotic behaviour. The variables 
a, p are defined in terms of x and Q 2 by 

a=[\n{x /x)\n{t/t )] 1 ' 2 , p = 

where t = ln(Q 2 /A 2 ). The fact that the data lie approximately on universal, 
horizontal bands is a demonstration of the validity of double asymptotic scal- 
ing. A slight breaking of the scaling behaviour is evident at large and small p 
in the lower plot, where, respectively, Q 2 becomes too small and x becomes too 
large for the approximation to be valid. Here, however, the full NLO DGLAP 
prediction (indicated by the curves) gives a very good description. 



ln(a; /x) 
lnft/to) 



(12) 
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Figure 6: HI F2 data compared with the 'double asymptotic scaling' prediction 1121. The 
curves are the full NLO best fit and have p = 2, 3, 4 in the a plot, and a = 1 , 1.2, 1.4 in the 

p plot. 



2.3 Beyond NLO DGLAP at small x 

Why does NLO DGLAP evolution work so well? To attempt to answer this, 
we consider two types of correction which could in principle be important 
in particular regions of x, Q 2 space. First, we note that the statement that 
'perturbative QCD evolution describes the data down to Q 2 ~ 1 GeV 2 ' must 
be qualified by the condition 'at small x\ At large x it has been clearly 
shown that higher-twist -/i.e. 1/Q 2 ) contributions to the structure functions 
are important. In Ref. u, a combined leading- twist /higher- twist analysis of 
structure function data, using the empirical form 

if* T W) = ^f T W)^, (13) 

gave C w 0.3 GeV 2 /(l — x) at large x and C « at small x. The factor 
(1 — x)^ 1 is not unexpected, since there are j times more operators which 
contribute to the jth moment of F 2 than to i 7 ^ • Therefore at large j 
the moments are expected to differ by a factor of j, corresponding to a factor 
(1 — a;) -1 difference-. in the contributions to the structure function. Although 
the analysis of Ref.ta was restricted to fixed-target data with x > 0.1 it would 
appear that the small-x HERA structure functions with Q 2 > 1.5 GeV 2 are 
also relatively free from higher-twist contributions. In fact a global fit to the 



10 



HERA (Q 2 > 1.5 GeV 2 ) F% data which includes a phenomenological higher- 
twist contribution of the form (a 2 /Q 2 )F2 LJ " > gives a value for a 2 consistent 
with zerollj. 

P(x,a s ) = £ nal [a/P n (x) 
xP n (x) = £ niSn . 1 A nm [ln(l/x)] m + xp n (x) 
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Figure 7: Pattern of \n(l/x) coefficients in the splitting functions at small x. 



A second important correction to NLO DLGAP evolution at small x 
comes from higher-order terms on the right-hand side of (||). As x — > 0, 
large L x = ln(l/x) logarithms appear in the splitting functions and spoil the 
convergence. The theory and phenomenology of these logarithms have been 
extensively discussed at this Workshopc3. In general one can show that there 
are at most n — 1 large logarithms in the nth order splitting function, i.e. 

xPW (x) = A n , n _i L™- 1 + 4,„-2ir 2 + ■ • • + A n , + ipW (x) , (14) 
where is finite in the limit x — » 0. The full splitting function is then 



xP(x, a s 



n=l 



V27T 



n-1 

E 

,m=0 



A n , m U x n + xp^(x) 



(15) 
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This expression holds for all four splitting functions, P gg , P qg , P gq and P qq . The 
pattern of the coefficients is illustrated in Fig. |7[ The horizontal rows represent 
the complete splitting function at a given order in perturbation theory, for 
example LO(n = 1), NLO(n = 2), etc. The leading diagonal with n = m + 1 
represents the set of leading L x contributions, and so on. 

Of particular phenomenological importance is the P qg function, since at 
small x <9-F2 /<9 In Q 2 ~ a s P qg ® g. To illustrate the importance of the higher- 
order contributions, it is convenient to take moments, j(u>) — J Q dx x u P(x), 
so that xP(x) ~ => 7(w) ~ r(n)/w™ The leading behaviour as u — > 

(equivalently x — > 0) of 7" is found to be Ell (with a s = 3a s /Tr) 

>f( U ) =n,f [(1 + .. .) + &. (^ + ---)+^ (^ + - 

' (16) 

where at each order in perturbation theory subleading terms down by one or 
more powers of to have been suppressed. In fact beyond NLO only the leading 
a s (a s /uj) n coefficients are known. In the limit a s L x 3> 1 they resum to give 
the characteristic BFKL 'perturbative pomeron' behaviourEil 

xP qg ~ x - iln2& ° . (17) 

This asymptotic limit is, however, not relevant for the HERA structure func- 
tion data. The interesting question is how important in practice are the con- 
tributions beyond NLO in ( |l6| ) . Certainly the leading logarithm part of these 
contributions appears to be large. This is illustrated in Fig. |^, which compares 
the exact leading- and next-to-leading-order contributions with the leading- 
logarithm part of the first five orders in perturbation theory. It would appear 
from Fig. || that the higher-order terms could be phenomenologically impor- 
tant. On the other hand, retaining only the leading terms at each order (as 
in the phenomenological analyses performed so far) could well overestimate 
their importance. This is evidently true of the LO and NLO contributions 
in the region 0.1 < u> < 1. Certainly, there is no evidence that the data re- 
quire such contributions. One can quantify this statement E3 by replacing 
\n(\/x) in the leading-logarithm expansion by \n.{xo/x)0(xQ — x), and regard- 
ing xq as a parameter to be determined by the data. Note that this replace- 
ment has the effect of artificially introducing sub-leading logarithms, since 

1.HERA F 2 
3). In each 



L% -► L% + nln(a;o)£™- 1 + • . . . Fig. | shows the X 2 of the fit to, 
data as a function of xq (in various resummation schemes, see 
scheme the minimum value of \ 2 is attained for .to — > 0, indicating that the 
data prefer fits without the higher-order leading ln(l/x) contributions. 
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Figure 8: The anomalous dimension ■y q9 (uj), showing the exact first- and second-order con- 
tributions, together with the leading-logarithm parts of the first five orders. 



24 'BFKL ' description of F 2 



In the double asymptotic limit discussed in Section 2/2, the evolution equation 
sums leading powers of [a s ln(l/a;) In Q 2 ] n generated by multigluon emission, 
and the distributions increase faster than any power of ln(l/:r) as x — ► 0, see 
( |ll| ) . The dominant region of phase space is where the gluons have strongly- 
ordered transverse momenta, Q 2 ^> k^ n ^> ... ^> k\ x . However such evolution 
does not include all the leading terms in the small- a; limit. It neglects those 
terms which contain the leading power of ln(l/cc) but which_are not accom- 
panied by the leading power of InQ 2 . The BFKL equationo, on the other 
hand, sums the leading ln(l/x) terms while retaining the full Q 2 dependence. 
The integration is taken over the full phase space of the gluons, not just 
the strongly-ordered part. The result is most conveniently established for the 
gluon distribution unintegrated over kr, 

rQ 2 

xg(x,Q 2 )= / dk 2 G(x,k^), (18) 
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Figure 9: x °f various resummed ln(l/x) fits to the HERA F2 data versus xq, from Ref.tfl. 

and is 

g{ x ,k 2 T ) ~ h{k 2 T ) x - x , (19) 

where h ~ (fc^) -1 ^ 2 at large fcy and A is the the maximum eigenvalue of the 
kernel K of the BFKL equation 



eg 

dL x 



d 2 k' T K(k T , k' T ) g(x, k'S) . (20) 



For fixed a s , A = 12 In 2a s /TT 2^.0.5. The prediction for F 2 is obtained by using 
the fcT-factorization theorem Ea 



F 2 (x,Q 2 ) = £^ J ^g^,k 2 T )F 2 (x',k 2 ,Q% 



(21) 



where F 2 denotes the quark-box contribution 717 — » qq for the scattering of a 
photon of virtuality Q 2 off a gluon with longitudinal momentum fraction x' 
and transverse momentum fcy- 

In recent years there have been several numerical analyses of the solution 
of the BFKL equation and of the corresponding predictions for F 2 , based on 
the above results (for a review seet3). The 'naive' BFKL prediction, 

F 2 (x,Q 2 ) = C(Q 2 )x- x + F? p , (22) 

where A « 0.5 and F^ F is a non-perturbative background contribution which is 
constant at small x, would appear to give too steep a rise in x. However there 
has been significant recent progress, reported at this Workshop c3, towards 
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KMS predictions for F 2 (x,Q 2 ) 
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Figure 10: Constrained CCFM description of F2 at small 



a 'unified' treatment which incorporates both DGLAP and BFKL dynamics, 
embodied in the so-called CCFM equation, together with appropriate kirie.- 
matic constraints. A reasonable description of the HERA data is obtaincdc3, 
see Fig. 10. Here the continuous curves are the CCFM predictions with and 
without the kinematic constraints. From a theoretical point of view, the at- 
traction of the BFKL approach is that it attempts to explain the shape of 
the small- x structure function from first principles. However, at present one 
cannot discriminate between the DGLAP and BFKL predictions on the basis 
of the x and Q 2 dependence of F2 alone. Furthermore, it is not clear how much 
truly 'BFKL' dynamics remains in the kinematically-constrained CCFM equa- 
tion applied to the HERA kinematic region. A key test would appear to be 
provided by more exclusive quantities, such as the average transverse hadronic 
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energy in deep inelastic events (which is predicted to be larger in the BEKL 
approach) and the cross section for producing forward, moderate Et jetsc3. 



2.5 Gluon determinations 

Deep inelastic structure functions directly measure the quark distribution func- 
tions. The precision is set by the experimental measurements, and is at the 
level of a few percent over a wide range in x, at least for the u and d quarks. 
The gluon is much less well determined: a variety of hard scattering processes 
provide measurements in particular x regions at different Q 2 scales, and mo- 
mentum conversation pins down the first moment to a few percent. 

Several new results on the gluon determination have been reported at 
this Workshop. An overview is given in Rcf. t3. We will concentrate here on 
those new results which are 'theoretically precise', i.e. are performed using 
cross sections calculated to next-to-leading order. This allows the gluon to 
be defined and extracted according to a particular factorization scheme (e.g. 
MS), as for the quark distributions. 

For x < 0.01, the Q 2 dependence of F2 is dominated by the gluon contri- 
bution to the right-hand side of the evolution equation, 



dF 2 (x,Q 2 ) a s (Q 2 )^ 2 f l dy 



The uncertainty in the gluon obtained from the HERA data in this way has 
steadily decreased in the last few years. As_an example, Fig. [ll] shows the 
gluon obtained from the HI NLO QCD fitE3. The shaded band represents 
the combined systematic and statistical error - approximately ±20% over the 
range in a; (10 -4 — 10~ 2 ) where there is sensitivity. It is interesting to compare 
the gluon of Fig. [ll] with that obtained from a global-fit. Fig. [[2] shows the 
gluons at 5 GeV 2 from the latest MRS(R) series of fitsE 2 !, together with several 
'old' gluons. At small x (~ 10 -4 ) there is good agreement with the HI gluon, 
and the spread of the new MRS gluons is similar to the uncertainty band in 
Fig. O. However for x > 10~ 2 the MRS gluon distributions are systematically 
smaller than the HI gluon. This is due to the extra constraints (principally 
from momentum conservation and from prompt photon production) imposed 
in the global fit. 

There has been considerable discussion at this Workshop about the new 
large Et jet data from CDF and DOlj. It is not yet clear whether the apparent 
excess of data over theory seen in the CDF data for Et ^ 200 GeV is a real 
effect, ancLifjSO whether it is due to harder parton distributions than previously 
cstimatedlld'Ell or to some 'new physics' contribution. What is clear is that the 
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Figure 11: The gluon distribution obtained by HI from a global NLO QCD fit 
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Figure 12: Gluon distributions at Q 2 
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jet data for Et ^ 200 GeV have considerable potential for constraining the 
gluon distribution (qg — > qg is the dominant subprocess) and a s . The CTEQ 
collaboration have already incorporated the CDF jet data into their global 
fits E3. As an example, Fig. |l^ shows a series of CTEQ gluon distributions 
obtained assuming different a s (M§) values in the range [0.105,0.122]. The 
effect of the jet data is seen in the difference between the new gluons and the 
CTEQ(3M) gluon in the x ~ 0.1 — 0.3 range. Notice also the correlation with 
the a s value. From such analyses, one may conclude that the uncertainty in 
the gluon | has now been reduced to O(±10%) over a broad range in x, rising 
to O(±20%) at small x ~ 10~ 4 . This represents considerable progress! 




o-l 1 1 1 1 1 1 — i — i — i — i — 1 

icr 4 icr 3 io~ 2 icr 1 .2 .3 .4 .5 .6 .7 x 



Figure 13: New CTEQ gluon distributions at 5 GeV 2 from fits incorporating CDF jet data, 13 !. 

The gluon distribution can also be determined at HERA from the '2 + 1' jet 
cross section, i.e. the cross section for producing an additional hard hadronic 
jet in addition to the current and beam remnant jets. Preliminary results had 
already been reported in Paris last year, but what is new this year is that the 
analysis can now be done consistently to next-to-leading order in perturbation 
theory. The key here is the_recent availability of NLO QCD Monte Carlo 
programs such as PRO JET E3 and MEP JET £3. At leading order the cross 
section is, schematically, 

<r(2 + 1) = fe ® ^(7*<? -» qg) + 9® °ii*g -> m)] ■ (24) 

Note that the subprocess cross sections are 0(a s ), and in fact this allows a 

c We refer here to the uncertainty at the starting scale, i.e. Qq ~ 4 — 5 GeV 2 . Perturbative 
evolution tends to make the distributions converge at higher Q 2 . 
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reasonably precise determination of the strong coupling, see Section § The 
parton distributions in fl24| ) are probed at x = XBj(l + Mj^/Q 2 ), where Mfj is 
the invariant mass of the final-state jet pair. In practice, the range covered at 
HERA is approximately 0.01 < x < 0.1. The importance of this is that it 'fills 
in' the gap between the small- a; scaling violatiop_.and large- ai_prompt photon 
methods for determining the gluon. Both the Hlo and ZEUSt£l collaborations 
have reported new results on the gluon distribution from er(2 + 1) at this 
Workshop. Fig. |lj shows the gluon distribution at 20 GeV 2 extracted by HI. 
The band corresponds to the combined systematic and statistical error. Note 
that there is good agreement with the standard gluons obtained from global 
fits. With a further reduction in the statistical and systematic errors, this 
method of determining the gluon will play an important role in the overall 
constraint picture. 
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Figure 14: HI 2 + 1 jet gluon distribution! 



Finally, we should mention also a new method of determining the gluon 
distribution from diffractive 3/tp electro- (and photo-) production: e + p — > 
e+J /ip+p. Fot^ufficiently high W 2 the scattering amplitude can be calculated 
perturbatively E3, reducing (at lowest order) to the scattering of the cc system 
off the proton via two gluon colour-singlet exchange. The lowest-order result 
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IS 



da 
~db 



(7*p^ J/ipp) 



T ee Ml^a s {Qy 
48aQ 6 



xg(x,Q ) 



Ml 



(25) 



where 



Q = (Q 2 + Ml)/ A , 



AQ /W 



(26) 



At higher orders, and when x is small, two gluon exchange is replaced by 
a generalized 'BFKL' gluon ladder, i.e. the unintegrated gluon distribution 
G(x, k 2 -) discussed in Section 2.4, but the basic structure remains the samecZI. 
Phenomenological studies-based on the above theoretical approach have re- 
cently been performed The W 2 dependence of the cross section is a 
direct measure of the shape of the gluon distribution, and the sensitivity is 
enhanced by ±he fact that it enters squared in the expression for the cross sec- 
tion. Fig. compares the predictions based on various gluon distributions 
with recent fixed-target and HERA data. Note that the theoretical calcula- 
tions do not yet include the complete 0(a s ) corrections to the cross sections, 
and cannot therefore distinguish gluons corresponding to different factorization 
schemes. In view of the potential of this method for constraining the gluon at 
small x, the calculation of the NLO corrections would be very worthwhile. 



T 




W (GcV) 



Figure 15: DifTractive J/ip photoproduction cross section data cormmred with theoretical 
predictions based on various gluon distributions EJ. 
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2.6 High Q 2 

Over most of the measured Q 2 range in Fig. El the deep inelastic cross section 
is dominated by the neutral current contribution, Fig. |l6|(a). With increased 




(a) (b) (c) 

Figure 16: (a) Neutral current, (b) charged current and (c) contact interaction contributions 
to electron— quark scattering at high Q 2 . 

statistics at high Q 2 ( > 10 4 GeV 2 ), HERA can also measure the structure of 
the proton as seen by the charged weak current, Fig. |l^(b). The corresponding 
deep inelastic cross section is 

dxdQ z 8snTcV(<3 + M^) 2 ^ 

(27) 

where A e is the helicity of the electron, Ui and refer to up- and down- 
type quarks respectively, and the V Ui d4 are the elements of the CKM matrix. 
In principle this allows a different combination of parton distributions to he 
measured. Some preliminary results have been reported at this Workshop El. 
The data clearly show the convergence of the charged and neutral current 
cross sections at very high Q 2 , da cc /dxdQ 2 ~ da NC /dxdQ 2 ~ Q -4 , see 
Fig. [l7| The fact that the high-Q 2 cross sections agree with the DGLAP-based 
extrapolations enables limits to be placed on possible contact interactions, as 
expected— on general grounds in composite models. The convention (see for 
examplel 41 ]) is to parametrize such interactions using the Lagrangian 

AC = ^ Ml,r , (28) 

q,l 

which leads to additional contributions to the cross section of order Q -2 A~ 2 
and A . The parameter A is a measure of the compositeness scale, and r\ 
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determines whether the interference with the Standard Model cross section 
is constructive or destructive. An example o£ how the data can be used to 
constrain A is shown in Fig. |l^ from ZEUS □. Here a particular (left-right) 
helicity structure has been assumed for the lepton-quark contact interaction. 
The actual lower limits on A depend slightly on this structure, and are typ- 
ically 0(1 TeV) and 0(2 TeV) for destructive and constructive interference 
respectively. Note that somewhat higher limits (3 — 4 TeV) are obtained from 
DrelLYan (qq — > l + l~) lepton-pair production at the Fermilab Tevatron col- 
liderc3. 



3 a s measurements 

Deep inelastic scattering has traditionally provided some of the most precise 
measurements of the strong coupling a s . There are three basic methods: 

• scaling violations of structure functions, 

• structure function sum rules, 

• jet cross sections. 
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Figure 18: Limits on lepton-quark contact interactirms from high-Q 2 neutral current DIS 

data, from ZEUSB. 

Some of the high-precision values |] obtained using these methods are listed in 
Table|^. The errors are between 5% and 10%, with the highest accuracy coming 
from the scale variation of the structure functions. More importantly, all the 
results agree within errors. These 'deep inelastic' a s measurements also have 
a significant weight in the overall world average value, which is currentlyEJ 

a s (M§) = 0.116 ±0.005. (29) 

The value of a s obtained from scaling violations has been discussed widely 
at this Workshop, especially in connection with new information coming from 
the HERA F% data. It is important to recall that the 'gold-plated' values in 
Table ^ come fijom the high-precision fixed-target data. In particular, in the 
analysis of Ref. which includes lower energy SLAC data and a phenomeno- 
logical contribution from higher- twists, the resolving power on a s comes from 
data with 0.3 < x < 0.5 (see Fig. |l|) and 10 < Q 2 < 100 GeV 2 . In this kine- 

^Throughout this section, a s refers to a s (AfJ,) in the MS renormalization scheme with 
n f = 5. 
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measurement 



experiment 



(Q) P QCD af s (Ml) 



slac/bcdmsi 
ccfrS 

CCFrII 

SLAC/EMC/SMcE 
HI, ZEUsEl 



dF l 2 N /d\nQ 2 
dF£ N /dlnQ 2 
tidx(F 3 » N + Fn 

o-7*p(2 + 1 jet) 



7.1 


NLO 


0.113 ±0.005 


5.0 


NLO 


0.111 ±0.006 


1.73 


NNLO 


n ins+ 006 
U.1U8_ _ 009 


1.58 


NNLO 


n i 99 +0.005 
U.1ZZ_ Q oog 


22.1 


NLO 


0.118 ±0.009 



Table 2: Some deep inelastic scattering a s measurements. 



matic region, higher-order perturbative and higher-twist corrections should be 
very small. 

In the last year it has become clear that the HERA structure function 
data prefer a slightly larger a—value than the fixed-target data at higher x, a 
result first pointed outin Ref.E3. For example, in the NLO analysis of the 1993 
data reported in Ref.B, the value a s LM 2 z ) = 0.120 ± 0.005(exp.) ± 0.009(th.) 
was obtained (for an update see Ref.ca), with the theory error dominated by 
the variation with respect to the renormalization and factorization scale. In 
the new MRS analysis 13, which includes the new 1994 HERA data, fits are 
performed with two a s values, 0.113 and 0.120. The latter value is motivated 
by (a) the preference of the LEP/SLC three-jeLand total Z hadronic width 
data for a larger a s value (see for example Ref. M), and (b) the fact that the 
Fermilab inclusive medium-_ET jet data is slightly steeper than the theoretical 
prediction using-jiartons with a s ~ 0.113 $3i3. Table |J shows the \ 2 values 
for the MRS fitsllj to the HI and ZEUS data with the two different a s values. 
Both experiments show a clear preference for the larger value. 



expt. 


x 2 (Ri) 

a s = 0.113 


x 2 (R 2 ) 

a s = 0.120 


HI (193 pts) 
ZEUS (210 pts) 


182 
391 


168 
362 



Table 3: \ 2 values for the new MRS(Ri,R 2 ) fits to the HERA F 2 datal 



What could be responsible for this apparent difference in the small- and 
large- a; a s values (if one ignores the fact that within the overall errors there 
is no significant disagreement!)? So far there has been no proper study of 
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Figure 19: Logarithmic Q 2 derivative of the F2 structure function measured off hydrogen 
and deuterium targets compared to NLO predictions, from Ref. 0. 



the dependence of the a s value on the way the charm contribution to F2 
is treated. This is almost certainly insignificant at large x, but could make 
a slight difference at small x where the charm contribution is expected to 
be a more significant fraction of the total structure function, and where the 
threshold Q 2 dependence of F$ could perhaps 'fake' a stronger DGLAP Q 2 
variation. Second, we have already seen that higher-order contributions to 
DGLAP evolution become more and more important as x — > 0. It could well 
be that the larger a s value at small x is an indication that the (as yet unknown) 
NNLO corrections are not negligible. A complete calculation of the third-order 
splitting functions would almost certainly settle this issue. 

New analyses of a s from the deep inelastic 2 + 1 jet cross section (Eq. (p4j)) 
by HI and ZEUS have also been reported at this Workshop It_appears 
that there may have been a slight problem with the previous analyses lj arising 
from a mismatch of the jet algorithm used in the NLO theoretical calculation 
and that employed in the experimental analyses. In general, jets are defined 
using a cluster algorithm which is applied to hadrons and partons in the ex- 
perimental and theoretical analyses respectively. One defines a dimensionless 
'metric' between two hadrons (say), = sy/Af 2 , and merges the hadrons 
into a single cluster if < y C ut, where y cu t is a fixed (small) number. Various 
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choices for (= 2EiEj(l — cos Oij) [JADE], = min fc^- [Durham/fcT], ■ ■ • ) 
and M 2 (= Q 2 ,W 2 , . . .) are possible, but in each case i/ie theoretical calcula- 
tion must be matched to the experimental definition — the NLO corrections are 
in general different for the different jet definitions. The new MEPJET theory 
Monte CarloEj allows an arbitrary algorithm to be used, and can therefore be 
tuned to the experimental analysis. Although there are still some questions to 
be answered (the effect of hadronization, the correlation of a s with the gluon 
distribution, the sensitivity to the jet cuts, etc.), it is already clear that this is 
a potentially powerful method for measuring the strong coupling at HERA in 
the future. It will be very interesting to see if the value thus obtained agrees 
with the scaling violation value. | 



4 Spin physics 

Polarized deep inelastic scattering experiments provide information about how 
the spin of the hadron is shared among its parton constituents. As for unpo- 
larized scattering, the cross section can be related to two structure functions, 
<7i and <?2- The former is analogous to F\, and in the parton model is given by 
a sum over polarized quark distribution functions: 

gi(x) = ±Y, e li A< i( x ) + A ^ • ( 3 °) 
i 

with Aq = q' — q+. In perturbative QCD, the distributions Aq and Ag ac- 
quire a scale dependence given by the polarized versions of the DGLAP equa- 
tions (ADGLAP). The first moments, J dxAqi and dxAg, represent the 
net spin carried by the various partons. The former can be related to axial- 
vector couplings measured in hyperon decay. In addition, the Bjorken sum 
rule (Jq dx(g^ — </") = glfl^i/ffylfl + 0(a s )]) provides a fundamental test of 
the theory. 

There are currently several important theoretical issues: 

• How much of the spin is carried by each type of parton; in particular, 
how much is carried by strange quarks and gluons? 

• Is the observed Q 2 dependence of the structure functions consistent with 
ADGLAP evolution? 

• How large is 172 and is it consistent with model calculations? 

e There is an analogy here with the two methods (total hadronic width and jet rates) for 
measuring a s at LEP/SLC. 
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• Can the theory predict the x — > behaviour of the polarized distribu- 
tions, as for unpolarized distributions (see Section pL4|)? 

Answers to all of these questions require precision input from experiments. The 
last few years have seen a rapid growth-.in the amount of data available. At 
this Workshop, the SMC collaborationEj have reported new measurements of 
gf, g2 and the semi-inclusive process fj,N — > NhX. The SLAC E-143 collabo- 
ration0 have reported new measurements of gi and the Q 2 dependence of g\' d . 
But perhaps the most significant development has been the first results from 
the HERMES experiment at DESYEE Fig. ||(] shows preliminary results on 
g™. The results are in agreement with published SLAC-E142 measurements Ej 
which cover approximately the same kinematic range. One of the advantages 
of the HERMES experiment is the ability to measure semi-inclusive processes, 
i.e. e^pTi e + h + X. The cross section for this is, schematically, 

Aa^, z) ~ &q(x)D Aq ^h(z) • (31) 

9 

The identification of the charge and flavour of hadrons in the final state there- 
fore gives a more powerful handle on the various Aq, in particular the separa- 
tion between valence and sea quarks. Quantitative results from HERMES are 
expected soon. 

New theoretical developments have also been reported at this Workshop. 
The major advance in the last year has been the calculation of the two-loop 
{0(a 2 )) polarized splitting functions^, 

AP ab (x,a s )^a s AP^\x)+a 2 s AP^(x) + ... . (32) 

This allows a consistent NLO ADGLAP analysis of the polarized structure 
function data. Already several groups iajzn_nerforrned global fits to determine 
the parton distributions Aqi and AgE§E§Ell. While the data can indeed be 
described by a simple set of starting distributions, there is not yet sufficient 
experimental information to determine the full quark flavour structure. With 
data on g\ and g™, only the valence Au and Ad quark distributions are well 
constrained. Further assumptions have to be made to extract the (Aqs) sea 
quark and gluon distributions (for example, SU(3) flavour symmetry for the 
sea, and Regge behaviour for the x — > behaviour of Aqg and Ag). As for 
the unpolarized distributions, weak information on the gluon distribution is 
obtained from the Q 2 dependence of the structure functions. Fig. ^l] shows 
anj-Gjxample of polarized distributions at Q\ = 4 GeV 2 obtained from a global 
fit e2J. The three sets of lines labelled (A), (B), and (C) correspond to fits with 
different assumptions about the form of Ag at large x. 
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First HERMES Results 



Spin Structure Function gj(x) of the Neutron: 
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Figure 20: The spin structure function g*(x) from HERMES E 



For the future, more data are urgently needed to further constrain the 
parton distributions. Improved precision on the Q 2 dependence of g\ at small 
x will help to determine Ag and the semi-inclusive measurements of and K~ 
will constrain Aqv and As respectively. The measurement of Ag at medium 
and large x presents a real challenge. As for the unpolarized distribution, what 
is needed is a process in which the gluon enters at lowest order. In this context, 
processes such as 



gtlptl 

Ti pT-l 
pl ipl i 



e + jets + 
J/?/> (or cc) 
(7, jets, . . 



X 



X 



X 



have been discussed at the Workshop, see for example Ref. An important 
point to remember is that for any such process, one must be sure that the 
collision energies, luminosities, cuts etc. are high enough that the unpolarized 
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gluon distribution can be reliably extracted from the corresponding unpolar- 
ized scattering process. 

Finally, the behaviour of g\ as x — > has been discussed T Jp_.tiLe context 
of a 'BFKL-type' resummation of large ln(l/x) logarithms Be3l1I. As for 
the unpolarized (singlet) structure function, the behaviour of (singlet) gf at 
small x is controlled by the i-channel exchange of a generalized gluon ladder. 
However there are important qualitative differences in the resulting behaviour. 
In the unpolarized case, both t-channel gluons are longitudinally polarized and 
the resulti ng b ehaviour is F 2 s (x, Q 2 ) ~ xg ~ (l/x) 0( " s ) ~ (1/x) 5 , as shown 
in Section 2.4. In the polarized case, this leading polarization configuration 
cancels, and the leading behaviour has one gluon transversely polarized. The 
predicted small- a; behaviour, obtained from resumming the terms proportional 
to [a s ln 2 (l/a;)]" to alLorders, is now the less singular gf(x, Q 2 ) ~ (1 / x)°(^°\ 
A detailed calculation^ gives a value A = 3.450Va s /(27r) w 1.0 for the power 
of 1/x. It will be interesting to see whether such behaviour can be measured 
experimentally. The same calculation gives the leading ln(l/x) contributions 
to the polarized splitting functions at all orders in perturbation theory, cf. 
Section 2^ . The impact ofpthese on the standard ADGLAP evolution at small 
x has been studied in Ref.E3. As for the unpolarized case, it is found that, in 
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practice, subleading logarithms are likely to play an important role. 



5 Photoproduction 

The high collision energies available at HERA allow photoproduction processes 
to be studied at very short distances. Hard scattering processes such as jet 
and charm production yield tests of perturbative QCD at next-to-leading order 
and provide information about the quark and gluon structure of the photon. 
Since the very first measurements of photoproduction cross sections at HERA 
in 1993, progress in this area has been steady. The new results of the past 
year are summarized in the report of the Photoproduction Working Group in 
these Proceedings. Here I will briefly mention two of these results which are 
particularly interesting. 

Large Et dijet photoproduction proceeds at leading order via the 'direct' 
scattering processes a^q — ► qg)®q p and 17(73 ~~ * QQ)^9 P j an d at higher orders 
through the 'resolved' processes a(qq — ► qq) ® q 1 <g> q p , cr(gq gq) (g> g 1 <g> q p , 
etc. The latter are particularly important at small Et- At leading order, the 
momentum fraction of the photon constituent can be reconstructed from the 
transverse energies and rapidities of the dijets, 

E T1 exp(-?7i) + E T2 exp(-r? 2 ) , on , 

x l = £gT ■ ( 33 ) 

Thus direct scattering events have x 7 « 1 and resolved scattering events have 
x 1 < 1. 

The angular distribution of the dijets probes the structure of the hard 
scattering process. In particular, small angle scattering is sensitive to the 
spin of the i-channel exchange particle. In this respect there is an important 
difference between the direct and resolved scattering cross sections, which are 
mediated by spin— ^ (quark) and spin— 1 (gluon) exchange respectively. Thus 

da dil _ 1 

— - g p <g> a(~fg ->■ qq) + . . . — ► : — 

dcostr 1 — I cos 0*1 

da IQS 1 

q 1 ®q p ® &{qq -+qq) + ... — > ; ^-73 , (34) 



dcosd* ' ' (l-|cos6l*|f 

where 9* is the scattering ang 



in the dijet centre-of-mass frame. In Fig. [22| , 
from the ZEUS collaboration E3, a cut in £ 7 (0.75) is used to separate the 
events into the two 'direct' and 'resolved' classes. The expected difference in 
the distribution as |cos0*| — > 1 between the 'direct' and 'resolved' samples 
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Figure 22: Dijet angular distribution, from Ref. 




is clearly seen, and in both cases there is good agreement with NLO QCD 
predictions. 

Another interesting result concerns the iirst quantitative measurements at 
HERA of large px photon photoproductionE^. At leading order, this proceeds 
via the Compton scattering of a photon off a quark: 

da ~ J2 4 1 P ® -> 7?) + • • • ■ (35) 

At higher orders there are also contributions from resolved processes such as 
qg — > jq. Fig. ^3] shows the distribution in x 1 for direct photon events from 
the ZEUS collaboration. A clear excess at large x 7 is seen, signalling the direct 
Compton scattering process. The interesting feature of the Compton process 
is that it is proportional to the fourth power of the quark charge, and therefore 
probes a different linear combination_pf quark distributions in the proton than 
the deep inelastic structure functionEJ. The photon angular distribution of the 
Compton process should again exhibit a characteristic (1 — | cos0*|) -1 form at 
small angles. 



6 Diffractive deep inelastic scattering 

The observation of 'diffractive' deep inelastic events with a large rapidity gapS 
has spawned a new field of interest at HERA in the last few years. At this 
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Figure 23: Distribution in x 1 of prompt photon events at HERA, from Ref. 



0. 



Workshop, as last year, there was much lively discussion about the measure- 
ment and interpretation of these events. 

The focus of attention so far is the diffractive structure function . This 
is obtained from a subclass of DIS events in which there is a large rapidity 
gap (HI, ZEUS) or from the observed excess of events with small Mx (the in- 
variant mass of hadrons in the detector) over the expected non-exponentially- 
suppressed non-diffractive contribution (ZEUS). The ZEUS collaboration have 
also identified energetic, small-angle protons in the Leading Proton Spectrom- 
eter (LPS), consistent with the interpretation that ep — > epX is the underlying 
process. 

The forward proton momentum can be parametrized in terms of the (small) 
momentum transfer t and the longitudinal momentum fraction 1 — x p ~ 1. 

One can then define a diffractive structure function F^* (xBj, Q 2 ', x p , t) or 
equivalently F° (4) (/?, Q 2 ; x p ,t) where (3 = x Bj /x p ~ Q 2 /(Q 2 + M%). If t 
is unmeasured and therefore integrated over, one can define a correspond- 
ing F® {A, Q 2 ; x p ) function. The crucial experimental observations are that 

i 7 "^ 4 ' is a leading twist structure function, indicating deep inelastic scattering 
off point-like objects, and that there is an approximate factorization property: 

F 2 D{3) ~ Xp - n T(P,Q 2 ) . (36) 

At the Paris Workshop last year, the values n = 1.19 ± 0.06(stat.) ± 0.07(sys.) 
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(Hl@) and n = 1.30 ± 0.08(stat.) +°; f 4 (sys.) (ZEUS@) were reported. This 
year, with new data, HI have measured n as a function of j3 and Q 2 , see Fig. g4j. 
A clear dependence on /3 (but not on Q 2 ) is observed. ZEUS have also reported 
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Figure 24: Measurement of the parameter n, where ~ Dj™, from HlEl. 

new values: n = 1.46±0.10 from defined by the Mx distribution0, and 

n = 1.28±0.17 from the LPS dataO. Note that these are not in contradiction, 
since the (/?) averages are different for the two samples. However, as last year, 
they do appear to be slightly larger than the HI values. 

The most attractive theoretical interpretation of these diffractive DIS events 
is that they correspond to the structure, function of the pomeron (V), rather 
than the proton (for a review see Ref. l3). In its simplest version, this model 
would predict 

f v (x p ,t)F?([3,Q 2 ) , 

i.e. a factorized structure function with n « 2a-p(0) — 1 w 1.16. This model is 
based on the notion ofjparton constituents in the pomeron' first proposed by 
Ingelman and SchleinO and supported by data from UA8. In such a model, a 
modest amount of factorization breaking could be accommodated by invoking 
a sum over Regge trajectories, each with a different intercept and structure 




F 2 DW (f3,Q 2 ;x pi t) = 
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function: 



F 2 D(4) (/3,Q 2 ;x p ,i) F^W^W- 1 f?(/?,Q 2 ) , 



(38) 



7? 



which would yield an effective n which depends on but is approximately 
independent of Q 2 . I— . 

The above approach has been put on a firmer theoretical footing in Ref.E-3, 
where the general concept of 'diffractive parton distributions' is introduced. In 
particular it is shown that an operator product definition exists, and that the 
diffractive distributions should satisfy DGLAP evolution equations. There has 
been a number of analyses based on DGLAP fits to diffractive structure func- 
tion data in the last year. At this Workshop the HI collaboration has presented 
a new measurement ofFi D ^ (i.e. integrated over both x-p and t) and a new 
NLO QCD fit, see Fig. |2^. The resulting parton distributions at Q 2 = 5 GeV 2 
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Figure 25: The diffractive structure function F^ ^ measured by HI I 

result of a NLO QCD fit. 



together with the 



are shown in Fig. g6|. Evidently the 'pomeron' (more precisely, the aggregate 
of the colour-singlet exchanges in Eq. (38)) is a predominantly gluonic object. 
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Although the structure function only directly measures the quark distributions, 
a very hard gluon distribution is needed to give an approximate scaling (in Q 2 ) 
behaviour at large /?, see Fig. ||[ 

From a purely phenomenological point of view, the above model has many 
attractive features. It enables predictions to be made, for example for the 
charm content of the diffractive .structure function Lj and the rate of jet pro- 
duction in diffractive DIS events L3. With a further assumption of universality 
of diffractive parton distributions, which lacks any theoretical proof at present 
however, one can make quantitative predictions for diffractive hard scattering 
processes in photon-hadron and hadron-hadron collisions EDO. 

Finally, we should mention a different approach to diffractive deep inelas- 
tic scattering based on perturbative QCD. In the microscopic colour dipole 
modelllj, a generalized 'BFKL' gluon ladder interacts with the qq Fock state 
of the virtual photon. In this model many features of can be predicted 
from QCD perturbation theory. The most striking features are the breaking 
of the Regge factorization of Eq. (|37|), and the qualitative differences between 
the longitudinal and transverse structure functions Lj: 

F? ~ 5 2 L ^ , FE ~ ±f (x p , . (39) 

Note that there is no x v — (J3, Q 2 ) factorization, and also no DGLAP evolution 
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as (3 — ► 1 at high Q 2 , where dominates. The x p dependence is determined 
by the x — > behaviour of the gluon distribution, and is generally steeper 
than the behaviour predicted by the simple (soft) pomeron model. For heavy 
quarks, a perturbatively hard scaLeJs set by m 2 , and infra-red safe predictions 
can be made from first principles I 79 " 39 ! Like the soft pomeron model, the pQCD 
model is so far in good agreement with the data. Once again the conclusion 
is obvious: more data, not only on F® {(3,Q 2 ]x pi t) but also on F® , charm 
and jet production, are urgently needed to discriminate between these different 
approaches. 

7 Conclusions 

At this Workshop we have seen significant progress in almost all aspects of 
deep inelastic scattering physics. Advances in experimental measurements are 
matched by advances in theoretical calculations, and our understanding of the 
short-distance structure of hadrons steadily improves as a result. However, 
each new Workshop inevitably produces new questions to be answered. Some 
of these have been discussed in this summary: for example, why does NLO 
DGLAP work so well even at small x and Q 2 values, can we push the un- 
certainty in the gluon distribution below the 10% level, is there a systematic 
difference between the deep inelastic a s values extracted from high and low 
x, how can the polarized gluon distribution be directly measured, and can the 
deep inelastic diffractive data really be understood in a simple 'pomeron par- 
ton model' or will the explicit factorization-breaking properties of the pQCD 
model be revealed? We look forward to some answers next year! 
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